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NUCLEOSYNTHESIS OF 26 AL IN ROTATING WOLF-RAYET STARS 
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ABSTRACT 



The 26 Al radionuclide can be detected through its de- 
cay emission line at 1.809 MeV, as was first observed 
by Mahoney et al. (1982). Since then, COMPTEL 
on board of the CGRO satellite, performed a sky 
survey in this energy range, and provided maps of 
the 26 Al distribution in the Galaxy. These results 
revealed that the main contributors to the synthesis 
of 26 Al are most likely the massive stars, which con- 
tribute through their winds (Wolf-Rayet stars) and 
through their supernova explosion. 
Comparison between these observations (in partic- 
ular observations in localized regions such as the 
Vela region and the Cygnus region) and the models 
available at that moment, showed however the need 
for improvements from both theoretical and observa- 
tional points of view, in order to improve our under- 
standing of the 26 Al galactic distribution as well as 
that of its synthesis. 

With the launch of the INTEGRAL satellite in Oc- 
tober 2002, the observational part will hopefully be 
improved, and the construction of better resolution 
maps at 1.809 MeV is one of the main aims of the 
mission. From a theoretical point of view, we need 
the most up-to-date predictions in order to be able 
to interpret the forthcoming data. 
In this paper, we address this latter part, and 
present new results for 26 Al production by rotating 
Wolf-Rayet stars and their contribution to the total 
amount observed in the Galaxy. 
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1. 26 AL SOURCES 

26 Al is a radioactive nuclide that can be produced 
by hydrostatic nucleosynthesis in H burning regions, 
by explosive nucleosynthesis and by spallation. Its 
half-life time in its ground state is of 7.2 x 10 5 yr, 



and it is thus a good tracer of recent nucleosynthctic 
events in the Galaxy. 

The maps obtained with COMPTEL (Diehl et al. 
1995, Knodlseder et al. 1999, Pliischke et al. 2001) 
allowed to pin down the main contributors to the dif- 
fuse emission observed mainly in the galactic plane 
: massive stars (Prantzos & Diehl 1996, Knodlseder 
1999). These objects end their lives with a supernova 
explosion, during which 26 Al can be synthesized and 
will eventually be expelled (Heger et al. 2003). The 
more massive (and shorter lived) ones will also con- 
tribute during quiescent evolutionary phases through 
their winds. 



2. WOLF-RAYET STARS 



Wolf-Rayet (hereafter WR) stars are the bare cores 
of massive O-type stars that lost their hydrogen 
rich envelope due to very strong stellar winds or by 
mass transfer in a close binary system. They exhibit 
at their surface the products of the internal nucle- 
soynthesis accompanying first the central H-burning 
(WN phase) and then central He-burning phase (WC 
and WO phases) . At solar metallicity, all stars more 
massive than about 30 M Q go through a WR phase, 
which typically lasts for a few 10 5 years in the case of 
non-rotating objects (Meynet & Maeder 2003; 2004), 
and massively contribute to the chemical enrichment 
of the interstellar medium (hereafter ISM) . 
26 Al is synthesized during the central H-burning 
phase in these stars, and is eventually ejected into 
the ISM as the star is peeled off by strong winds dur- 
ing the WR phase. The amount of 26 Al ejected in 
the ISM through winds is very sensitive to mass loss 
rates, metallicity, rotation and initial mass. Meynet 
et al. (1997) presented a first study of the effects 
of initial mass and metallicity, as well as the effects 
of different possible prescriptions for convection. At 
that time, they could already show that the amount 
of 26 Al ejected increases with metallicity and ini- 
tial mass. In the grid of models used by Meynet et 
al. (1997), artificially twice enhanced mass loss rates 
were used, which appeared at that time to better fit 
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the observed properties of the WR stars in a series 
of environments (Maeder & Meynet 1994). 

In the following, we present the results obtained for 
a new grid of models including rotation and updated 
mass loss prescriptions for the different evolutionary 
phases that massive stars are likely to go through. 
Part of these results were discussed in Vuissoz et al. 
(2004). 



3. ROTATION INDUCED MIXING 

Massive stars are known to have high equatorial ve- 
locities (<v> = 200 - 250 km.s -1 ) while on main 
sequence. Such a rotation affects the evolution of 
stars in several ways. It first has a direct incidence 
on its structure through the centrifugal acceleration. 
It also affects the mass loss in terms of geometry 
and quantity. Finally it triggers instabilities in the 
radiative interior of the star, which allow transport 
of angular momentum and chemical species. 

Transport of angular momentum and chemical 
species is ensured by the meridional circulation, 
which settles in to balance the thermal desequi- 
librium induced by the departure from spherical 
symetry in rotating objects, and by the shear insta- 
bility associated with the turbulent layers in regions 
with high angular velocity gradients (Zahn 1992; 
Maeder & Zahn 1998). 

We chose an initial velocity v = 300 km.s -1 for 
all masses and metallicities, taking into account that 
this value is compatible with the mean averaged ob- 
served velocity of O-type stars at solar metallicity. 
Considering this, the rotating models used here well 
reproduce a series of observed features for O-type 
and WR stars, that non-rotating models fail to fit 
(Meynet & Maeder 2004). 



4. RESULTS FOR 26 AL 

The effects of rotation on 26 Al production in WR 
stars are consequences of both the effects on the mass 
loss and of the rotation induced mixing. The main 
results for what concerns 26 Al are the following : 

• At given metallicity, rotation allows stars with 
lower initial masses to enter the WR phase. 
When we take stellar rotation into account 
(which is an observational fact), the number of 
stars contributing to the 26 Al enrichment of the 
ISM through stellar winds is thus increased. The 
total contribution of WR stars to the galactic 
content of 26 Al is larger when using rotating 
models than when using non-rotating ones, as 
can be seen in Fig. 3, and as will be discussed 
in the next section. 

• At given mass and metallicity, rotation results 
in an entrance into the WR phase at an earlier 
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Figure 1. Evolution as a function of time of the cu- 
mulated mass of 26 Al ejected by WR stellar winds 
at twice the solar metallicity. Continuous line show 
the case of rotating models, dotted lines the case of 
non-rotating models. The short-dashed lines and the 
long-dashed lines show the evolution as a function of 
time of the cumulated mass of "live" 26 Al , i.e. of 
still non-decayed 26 Al , for rotating and non-rotating 
cases respectively.. 
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Figure 2. Same as Fig. 1 for solar metallicity mod- 
els. 
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evolutionary stage (not necessarily at an ear- 
lier time). For instance the solar metallicity 60 
M star model enters the WR phase during the 
core He-burning phase when it is non-rotating 
and during the end of the core H-burning phase 
when it is rotating. In this last case the duration 
of the WR phase is longer so that the cumula- 
tive time interval during which 26 Al is expelled 
is larger for rotating stars. This leads to an en- 
hancement of the 26 Al yields. This is well illus- 
trated in Figs. 1 and 2 where the evolution of 
the cumulated mass of 26 Al is shown as a func- 
tion of time for both rotating and non-rotating 
models. 

• The higher the metallicity, the larger the ejected 
mass of 26 Al (see Fig. 1 and Fig. 2). 

This is due to the metallicity dependence of 
mass loss rates and to the higher 25 Mg mass 
fraction available for 26 Al synthesis at higher 
metallicity. 

• Rotation allows surface enrichment of 26 Al prior 
to the WR phase. This is a direct consequence 
of the combined effects of meridional circulation 
and turbulence which allow 26 Al to diffuse out- 
wards from the convective core, where it is pro- 
duced by proton capture on 25 Mg, up to the 
radiative envelope, and to appear earlier at the 
surface of the star. 



5. CONTRIBUTION OF WR STARS TO THE 
26 AL GALACTIC CONTENT 

One can evaluate the contribution of WR stars to 
the total galactic mass of 26 Al in the following way: 

rR /-120 
= / 2nrdra(r) / <t>(M)Y 2 ^ (Mi,r, Z)dM 

JO VM LWR (r,Z) 

if T 26 = 1.44 Myr » 1 Myr, (1) 

R is the galactic radius taken equal to 15 kpc. 
<j(r) is the surface density distribution of atomic and 
molecular hydrogen according to Scoville & Sanders 
(1987). This quantity is used to mimic the present 
star formation rate in the Galaxy, provided that the 
molecular clouds distribution is considered a good 
tracer of the star forming regions. 
Mlwr is the minimum mass for a star to enter the 
WR phase at a given metallicity. Its dependence on 
initial velocity and initial metallicity is explicitely 
taken into account. 

$(M) oc m~( 1+a: ) is the Initial Mass Function, x 
defining the slope of the IMF. Here, two different 
prescriptions were used, in order to show the impact 
of this physical ingredient on the global production 
of 26 Al by WR stars. 

Y$ {Mi, r, Z) is the yield of 26 Al (excluding the su- 
pernova contribution) for a given mass Mj at a given 
metallicity Z. 

The expression in Eqn. 1 is normalized so as to ac- 
count for a rate of 3 SNae per century in the Galaxy 



(e.g. Capcllaro et al. 1999). 

Furthermore, we considered a metallicity gradient 
along the galactic radius such as Z = 0.04 for r < 
4 kpc, and d log(Z)/dr = -0.07 dex.kpc -1 beyond, 
with the condition that Z = 0.02 at r = 8.5 kpc (so- 
lar neighbourhood). 

In Fig. 3, we can see that the radial surface density 
profile of 26 Al has a similar shape whether we use 
rotating or non-rotating models. In particular, in 
both cases there is a peak between 4 kpc and 6 kpc 
associated with the existence of the molecular clouds 
ring at 5 kpc, which is accounted for in the expression 
of the surface density distribution cr(r). 



• The rotating models lead to a larger contribu- 
tion of WR stars to the 26 Al surface density pro- 
file in the inner 8 kpc ring. There, Z > Z©, and 
as seen in Fig. 1 and Fig. 2, the combined effects 
of metallicity and rotation lead to a substantial 
enhancement of the yields. 

• The total mass of 26 Al originating from WR 
stars is increased by a factor 2.5 to 2.9 when 
rotation is taken into account. The values ob- 
tained are also quite sensitive to the IMF slope, 
the Salpeter IMF (x = 1.35) favouring a larger 
contribution, even though the Scalo IMF (x = 
1.7) is considered better suited for massive stars. 

• For an IMF slope of 1 .7, we obtain a total ejected 
mass of 26 Al by WR stars of 1.27 M Q in case 
of rotating models, which represents 42% to 63 
% of the total observed mass, that is estimated 
to range between 2 and 3 M©. 



According to the recent observations by RHESSI 
(a 7-ray observatory for the Sun which was able 
to measure the diffuse emission of both 26 Al and 
60 Fe), the new estimate derived for the line ratio 
60 Fe/ 26 Al ~ 0.1 (see Smith's contribution in these 
proceedings) is consistent with the value predicted by 
Timmes et al. (1995), from computations including 
the sole contribution of type II supernovae (SNII). 
This would designate SNII as the dominant contrib- 
utors to the 26 Al galactic content. However, the most 
recent computations of nucleosynthesis from SNII in- 
dicate that the line ratio 60 Fe/ 26 Al should be larger 
than about 0.4. If this is the CctSG, clS pointed out 
by Prantzos (2004; see also his contribution in these 
proceedings), 26 Al should, at least up to a 50 % level, 
originate from sources that produce 26 Al but do not 
produce 60 Fe, as is the case for the WR stars. Our 
results show that rotating stellar models of WR stars 
could actually be responsible for such an amount of 
26 Al . Let us stress however, that this result has to 
be taken with caution, considering the 2.6 sigma level 
of the RHESSI detection of 60 Fe and the large uncer- 
tainties existing on the physical parameters entering 
the derivation of such a quantity (see also Palacios 
et al., 2004). 
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Figure 3. Theoretical profiles of the radial surface density of 26 Al in the Galaxy (hatched histograms) compared 
with the observational profile obtained after COMPTEL data analysis (Knodlseder 1997) (black line). The total 
26 Al mass ejected by WR winds in the whole Galaxy during the last million years (in Mq) is indicated on each 
plot. Upper panels: Profiles obtained for rotating models v- m i = 300 km.s , for a "Scalo" IMF slope x = 1.7 
(on the left), and a "Salpeter" IMF slope x — 1.35 (on the right). Lower panels: Same as upper panels, but 
using non-rotating models. 



6. CONCLUSIONS 



We have presented new results concerning 26 Al pro- 
duction by WR stars from a new grid of models in- 
cluding rotation and updated physics, in particular 
the most recent prescriptions for mass loss rates. 
Taking rotation into account globally leads to an en- 
hancement of wind ejected mass of 26 Al by very mas- 
sive stars. 

Convolved with appropriate IMF and star formation 
rate indicator, these yields lead to a total mass of 
about 1.3 Mq of 26 Al originating from WR stellar 
winds. This value is in agreement with the conclu- 
sions drawn by Prantzos (2004) from recent measure- 
ments of the line flux ratio 60 Fe/ 26 Al. 
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